Search for: All records

ABSTRACT Semi-analytic modelling furnishes an efficient avenue for characterizing dark matter haloes associated with satellites of Milky Way-like systems, as it easily accounts for uncertainties arising from halo-to-halo variance, the orbital disruption of satellites, baryonic feedback, and the stellar-to-halo mass (SMHM) relation. We use the SatGen semi-analytic satellite generator, which incorporates both empirical models of the galaxy–halo connection as well as analytic prescriptions for the orbital evolution of these satellites after accretion onto a host to create large samples of Milky Way-like systems and their satellites. By selecting satellites in the sample that match observed properties of a particular dwarf galaxy, we can infer arbitrary properties of the satellite galaxy within the cold dark matter paradigm. For the Milky Way’s classical dwarfs, we provide inferred values (with associated uncertainties) for the maximum circular velocity $$v_\text{max}$$ and the radius $$r_\text{max}$$ at which it occurs, varying over two choices of baryonic feedback model and two prescriptions for the SMHM relation. While simple empirical scaling relations can recover the median inferred value for $$v_\text{max}$$ and $$r_\text{max}$$, this approach provides realistic correlated uncertainties and aids interpretability. We also demonstrate how the internal properties of a satellite’s dark matter profile correlate with its orbit, and we show that it is difficult to reproduce observations of the Fornax dwarf without strong baryonic feedback. The technique developed in this work is flexible in its application of observational data and can leverage arbitrary information about the satellite galaxies to make inferences about their dark matter haloes and population statistics. 

Abstract When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature.To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m= 50 cm²/g).When simulating these halos, we also vary the massresolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses.The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10⁶particles per halo), we find that variations of order 10% in collapse time are still present.The results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context. 

Abstract We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (f_diss= 0.75). In this paper, we examine the properties of dwarf galaxies withM_*∼ 10⁵–10⁹M_⊙in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm²g⁻¹. On the contrary, models with (σ/m) = 10 cm²g⁻¹produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm²g⁻¹are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm²g⁻¹,f_diss= 0.75 are in tension in massive dwarfs (M_halo∼ 10¹¹M_⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures. 

Abstract The properties of globular clusters (GCs) contain valuable information of their host galaxies and dark-matter halos. In the remarkable example of ultra-diffuse galaxy, NGC5846-UDG1, the GC population exhibits strong radial mass segregation, indicative of dynamical-friction-driven orbital decay, which opens the possibility of using imaging data alone to constrain the dark-matter content of the galaxy. To explore this possibility, we develop a semianalytical model of GC evolution, which starts from the initial mass, structural, and spatial distributions of the GC progenitors, and follows the effects of dynamical friction, tidal evolution, and two-body relaxation. Using Markov Chain Monte Carlo, we forward-model the GCs in a UDG1-like potential to match the observed GC statistics, and to constrain the profile of the host halo and the origin of the GCs. We find that, with the assumptions of zero mass segregation when the star clusters were born, UDG1 is relatively dark-matter-poor compared to what is expected from stellar-to-halo–mass relations, and its halo concentration is lower than the cosmological average, irrespective of having a cuspy or a cored profile. Its GC population has an initial spatial distribution more extended than the smooth stellar distribution. We discuss the results in the context of scaling laws of galaxy–halo connections, and warn against naively using the GC-abundance–halo–mass relation to infer the halo mass of ultra-diffuse galaxies. Our model is generally applicable to GC-rich dwarf galaxies, and is publicly available. 

Abstract Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. BelowM_⋆∼ 10⁹M_⊙, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way (MW). Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the largest sample of satellite galaxies in the Local Volume (D≲ 12 Mpc) to date. We use 250 confirmed and 71 candidate dwarf satellites around 27 MW-like hosts from the Exploration of Local VolumE Satellites (ELVES) Survey and use the semianalyticalSatGenmodel for predicting the population of dark matter subhalos expected in the same volume. Through a Bayesian model comparison of the observed and the forward-modeled satellite stellar mass functions (SSMFs), we infer the satellite stellar-to-halo mass relation. We find that the observed SSMF is best reproduced when subhalos at the low-mass end are populated by a relation of the form $M_{\star }\propto M_{\mathrm{peak}}^{\alpha }$, with a moderate slope of ${\alpha }_{\mathrm{const}}=2.10\pm 0.01$and a low scatter, constant as a function of the peak halo mass, of ${\sigma }_{\mathrm{const}}={0.06}_{-0.05}^{+0.07}$. A model with a steeper slope (α_grow= 2.39 ± 0.06) and a scatter that grows with decreasingM_peakis also consistent with the observed SSMF but is not required. Our new model for the satellite–subhalo connection, based on hundreds of Local Volume satellite galaxies, is in line with what was previously derived using only MW satellites. 

Abstract Isolated dwarf galaxies usually exhibit robust star formation but satellite dwarf galaxies are often devoid of young stars, even in Milky Way–mass groups. Dwarf galaxies thus offer an important laboratory of the environmental processes that cease star formation. We explore the balance of quiescent and star-forming galaxies (quenched fractions) for a sample of ∼400 satellite galaxies around 30 Local Volume hosts from the Exploration of Local VolumE Satellites (ELVES) Survey. We present quenched fractions as a function of satellite stellar mass, projected radius, and host halo mass, to conclude that overall, the quenched fractions are similar to the Milky Way, dropping below 50% at satelliteM_*≈ 10⁸M_⊙. We may see hints that quenching is less efficient at larger radii. Through comparison with the semianalytic modeling codeSatGen, we are also able to infer average quenching times as a function of satellite mass in host halo-mass bins. There is a gradual increase in quenching time with satellite stellar mass rather than the abrupt change from rapid to slow quenching that has been inferred for the Milky Way. We also generally infer longer average quenching times than recent hydrodynamical simulations. Our results are consistent with models that suggest a wide range of quenching times are possible via ram pressure stripping, depending on the clumpiness of the circumgalactic medium, the orbits of the satellites, and the degree of earlier preprocessing. 

ABSTRACT We combine the isothermal Jeans model and the model of adiabatic halo contraction into a semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. The model agrees well with cosmological SIDM simulations over the entire core-forming stage up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable feature in the context of the structural diversity of bright dwarfs. The fast speed of the method facilitates analyses that would be challenging for numerical simulations – notably, we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, Mb/Mvir, and galaxy compactness, r1/2/Rvir; we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and we delineate the regime of core-collapse in the Mb/Mvir − r1/2/Rvir space, for a given cross section and concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible improvements of the fluid model. We have made our model publicly available at https://github.com/JiangFangzhou/SIDM. 

ABSTRACT Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass haloes, accelerated if haloes are subhaloes of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of $${\lesssim}10^8 \, \mathrm{M}_\odot$$ subhaloes in lens-mass systems ($${\sim}10^{13}\, \mathrm{M}_\odot$$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhaloes with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhaloes (≲1/1000 host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhaloes to core-collapse, finding that it is nearly impossible to drive core-collapse in subhaloes in SIDM models with constant cross-sections. Any discovery of ultracompact dark substructures with future substructure lensing observations favours additional degrees of freedom, such as velocity-dependence, in the cross-section.